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We study the behavior of a semi-infinite monolayer, which is placed initially on a half of an infi- 
nite in both directions, ideal crystalline surface, and then evolves in time due to random motion of 
the monolayer particles. Particles dynamics is modeled as the Kawasaki particle-vacancy exchange 
process in the presence of long-range attractive particle-particle interactions. In terms of an an- 
alytically solvable mean- field- type approximation we calculate the mean displacement X(t) of the 
monolayer edge and discuss the conditions under which a monolayer spreads (X(t) > 0), partially 
wets (X(i) — 0) or dewets from the solid surface (X(t) < 0). 
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Dynamics and static properties of thin liquid films on 
solid surfaces have been studied for many years result- 
ing in a seemingly good understanding of the problem 
jl],^| . However, with the advent of new experimental 
techniques, capable of studying properties of molecularly 
thin (MT) films, it has become clear that the developed 
theoretical concepts apply only to sufficiently thick films; 
for MT films significant departures from the standard be- 
havior have been observed In particular, several 
remarkable features have been revealed by ellipsometric 
studies of the MT precursor films, i.e. films emitted by 
(sessile) liquid drops placed on solid substrates Q: 
First, such films have been detected even in the case of 
non- wetting drops. This implies that physical conditions 
at which such a MT film appears may be different of the 
ones corresponding to the wetting/dewetting transition 
at macroscopic scales. Next, precursors do not spread 
at a constant rate; the mean displacement of the film's 
edge grows with time t only in proportion to \/t. Lastly, 
" fine structure" of the MT precursors may be very differ- 
ent; in some cases the film's density shows a pronounced 
variation with the distance from the macroscopic drop, 
which reveals the surface-gas-like, rather than the liquid- 
like behavior. In other systems, the films are dense and 
compact. Even more striking, on the intermediate-energy 
substrates surprising "terraced" patterns appear, formed 
by several superimposed MT precursors each spreading 
at the yt-rate on top of lower layers. 

Meanwhile, several attempts have been made to ex- 
plain why do the MT precursor films spread at the yt- 
rate. Ref.5 proposed a "stratified droplet" model, in 
which a sessile drop is regarded as a succession of hor- 
izontal layers, each layer being a two-dimensional, in- 
compressible Navier-Stokes liquid. This model suggests 
that the yi-l&w results from the competition between the 
liquid-solid attractions, which represent the driving force 
of spreading, and viscous- type frictional forces, which 
control particles dynamics on the solid surface. We note 
parenthetically that similar ideas have been used to de- 



scribe dynamics of the reverse processes - dewetting of 
a monolayer B and squeezing of a MT lubricating film 
out of a gap between two solids @] , for which it has been 
also predicted that the radii of the dewetted areas grow 
at the Vt-rate. Next, Ref.8 described droplet spreading 
in terms of the Langevin dynamics of a non- volatile fluid 
edge, modeled by horizontal solid-on-solid-model strings. 
Such an approach has reproduced the "terraced" pro- 
files; the Vt-law was found, however, only as a transient 
regime. Lastly, in a microscopic approach of Refs.9 the 
MT film was considered as a lattice gas of interacting 
particles connected to a reservoir (droplet). Here, the 
-\/t-law was obtained for both precursors of the sessile 
drops and creeping films in the capillary rise geometries; 
it was claimed that such a behavior is controlled by mi- 
gration of voids from the advancing edge of the film to 
the reservoir. 

Despite reasonably good explanation of the dynamical 
behavior, provided by Refs.5,8 and 9, several fundamen- 
tal questions still remain largely unanswered. In partic- 
ular, the dependence of the prefactor in the y/t-law on 
the temperature, keT, and on the parameters of the in- 
teraction potentials has not been elucidated so far. As 
a matter of fact, the model of Ref.5 discards the effects 
of the drop's surface tension and/or of the monolayer 
edge tension j e on spreading kinetics. In consequence, 
Ref.5 predicts that "terraced" spreading appears as soon 
as any kind of attractive liquid-solid interactions (LSI) is 
present, which contradicts apparently to the experience 
Q . Contrary to Ref.5, the models of Refs.8 and 9 take 7 e 
into account and show that the film may actually appear 
only if the strength of the LSI exceeds certain threshold 
value. However, a common subtle point of both Refs.8 
and 9 is that the prefactor in the Vt-law is expressed in 
terms of several parameters, which are assumed to be in- 
dependent of the dynamics; in Refs.9, for instance, these 
are the particle density in the reservoir and 7 e . On the 
other hand, 7 e originates from attractive liquid-liquid in- 
teractions (LLI) and thus depends on the density profile 
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in the film. The latter is itself dependent on the spread- 
ing rate and hence, on j e . Therefore, the calculation of j e 
and, consequently, of the prefactor in the \/t-law requires 
solution of essentially non-linear dynamical problem in 
which attractive LLI are taken into account explicitly. 

In this Rapid Communication we study analytically 
the behavior of a liquid monolayer, which occupies ini- 
tially a bounded, macroscopically large area of the solid 
surface, and then evolves in time due to random motion of 
the monolayer particles. Particles' dynamics is modeled 
as the Kawasaki-type particle-vacancy exchange process 
in the presence of short-range repulsive (hard-core) and 
weak long-range attractive particle-particle interactions. 
Here we consider a simple case when the initially occu- 
pied region is the half-plane — oo < X < 0, Fig.l, and 
calculate the mean displacement X(t) of the monolayer 
edge. We note that our results apply, as well, to the inter- 
mediate time behavior in several other two-phase geome- 
tries. Particularly, the initially dewetted region can be a 
hole of radius R, nucleated in a homogeneous monolayer, 
or the monolayer can occupy a circular region of radius 
R, which situation appears at the late stages of sessile 
drops spreading B, For such geometries, our results de- 
scribe the kinetics on time scales such that X(t) <C R, 
in which regime the precise form of the phase-separating 
boundary is not important (see, e.g. Refs.9). 

To determine the time evolution of X(t), we develop 
a mean-field-like, self-consistent approach, in which the 
non-linear coupling between the density distribution in 
the spreading film and the edge tension 7 e is taken into 
account explicitly. Within this approach we recover the 
result of Refs.9, i.e. the law X(t) — AyTJot, in which 
Do is the bare diffusion coefficient describing dynamics 
of an isolated particle on the solid surface. Here, how- 
ever, we define the prefactor A explicitly as a function of 
keT and of the interaction parameters. We show that A 
can be positive or negative, which means that the mono- 
layer can spread, partially wet or dewet from the solid 
surface, and determine the temperature T w / dw of the 
wetting/dewetting transition in the monolayer regime. 
Moreover, we find that spreading of the monolayer can 
proceed quite differently at different T, which agrees with 
the experimentally observed behavior [Q. When T > Tb, 
where Tb is also found explicitly, we have that "f e = 0, 
A ~ -y/ln(t) as t — > oo , and the density in the monolayer 
varies strongly with the distance from the edge. We re- 
mark that this finding contradicts to Refs.9, which sug- 
gest that such a "surface-gas" -like spreading may take 
place only in the absence of attractive LLI. For lower T, 
such that T w / dw < T < Tb, we find that the density 
variation is less pronounced and both A and 7 e are posi- 
tive and constant, which signifies that in this T-range the 
monolayer spreads as a " liquid" . Lastly, for T < T w /^w 
the monolayer dewets from the surface. 

The model to be studied here is defined as follows: 
(a) The particles experience two types of interactions - 
the LLI and the LSI. The LSI create effectively a lattice 
of potential wells (with the coordination number z and 



the interwell distance a) , such that the particles reside in 
the local minima of these wells. We assume that the LSI 
correspond to the limit of the so-called intermediate lo- 
calized adsorption, which is appropriate for many adsor- 
bates and persists over a wide T-range §]. In this limit 
the particles are neither completely fixed in the wells, 
nor completely mobile: The LSI wells are deep with re- 
spect to desorption (desorption barrier Ud keT) so 
that only a monolayer can exist, but have much lower 
barrier Ui against the movement across the surface. Fur- 
ther on, we suppose that the LLI are two-body, cen- 
tral and additive; the LLI potential U(r) is a hard-core 
at distance r = <x, which means that each well can be 
occupied by one particle at most, and is attractive for 
r > cr, U(r) = -U (T)(er/r) n , n > 2. The parameter 
Uo(T) <C Ud, which is the case for many realistic situa- 
tions J2| and which implies that the LLI incur only small 
local perturbations to the array of the LSI wells. 

(b) Occupation of the well with radius-vector r at time 
t is described by the variable 7j(r;t), which can assume 
two values, and 1; its realization average value, i.e. the 
local density, is denoted as p(r;t) = rj(r;t). The initial 
configuration of the monolayer is depicted in Fig.l, i.e. 
the monolayer particles are placed at random positions 
and at a fixed coverage p < 1 (number of occupied wells 
as a fraction of the total number of wells per unit area) 
in the wells of the half-plane — 00 < X < 0. 

(c) The particles motion is activated by chaotic vibra- 
tions of solid atoms and proceeds by rare events of hops 
between the local minima of adjacent wells. In absence 
of the LLI, one may estimate the diffusion coefficient for 
such a motion to be Do ~ a->er 2 /z, where u is the fre- 
quency of hops. Now, the LLI couple the dynamics of 
any given particle to the motions of all other particles: 
first, hard-core repulsion prevents multiple occupancy of 
any potential well; second, on escaping from the well with 
radius-vector r, any given particle follows preferentially 
the local gradient of the LLI potential landscape U(r;t): 

U(r,t) = -Uo(T) CT »^^^, (1) 

where the summation extends over the entire lattice, ex- 
cluding r" = r. To take the LLI into account, we stip- 
ulate that for any given particle leaving the well with 
radius- vector r at time moment t, the choice of the jump 
direction is random and governed by the position- and 
time-dependent probabilities |k|: 

p(r\r' z ) = Z- 1 exp(|[U(r;t)-U(r; ; t)]), (2) 

where r' z is the radius-vector of one of z wells, adjacent 
to the well at r, (3 = I/keT, and Z denotes the normal- 
ization factor. When the jump direction is chosen, the 
particle attempts to hop into the target well. The hop is 
fulfilled if this well is unoccupied; otherwise, the particle 
attempting to hop is repelled back to the well at r. 
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We proceed further on by assuming local equilibrium 
(see, e.g. Refs.10), which implies that occupations of dif- 
ferent wells factorize and thus allows for the description 
in terms of local densities, p(r;t). For these, we find 



~ §- t P(r;t) = - p(r;t) J2 F ,p(r\r> z ) (1 - p(r' z ,t)) 
+ (1 - p(r;t)) £ r 'PKk)p(rL;t), 
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in which p(r|r^,) are determined by Eqs.(2) and (1) with 
n(r; t) replaced by p(r; t). Eq.(3) has to be solved subject 
to the initial condition that p(r; t = 0) = for X > 
and p(r; t = 0) = p for X < 0. 

Let us analyse now, on the basis of Eq.(3), the time 
evolution of the mean displacement X(t) of the mono- 
layer edge. To do this, we follow Refs.9 supposing that 
for the long-ranged, but rapidly vanishing LLI, a hop of 
any particle which is not directly at the edge, does not 
change the energy in Eq.(l). This means, in virtue of 
Eq.(2), that for such particles all hopping directions are 
equally probable and their migration on the surface is 
constrained by the hard-core interactions only. The par- 
ticles being at the edge, however, are effectively attracted 
by the "bulk" monolayer, which results in asymmetric 
hopping probabilities (see Refs.9 for more details). We 
note that such an approximation is actually a translation 
onto the molecular level of standard descriptions of the 
liquid front dynamics as a competition between surface 
tension and internal pressure |l],|2| . 

Further on, since we are interested in the behavior 
of the mean displacement of the edge, it is justified to 
neglect the fluctuations around X(t) along the Y-axis. 
Consequently, we suppose that the edge position docs 
not depend on Y, Fig.l, which makes the system effec- 
tively one-dimensional and the original two-dimensional 
geometry enters only through the 2D diffusion coefficient 
D and 2D edge tension, (see Eq.(9)). 

Now, the one-dimensional model we have to study con- 
sists of a ID hard-core lattice gas, which is put initially 
into a "shock" configuration and then evolves in time by 
particles attempting to hop to the nearest unoccupied 
sites. The "shock" configuration means that all particles 
are placed at random with a fixed mean density p at the 
sites — co < X < of an infinite in both directions ID 
lattice. All particles, except the rightmost one (RMP) 
which defines position of the edge, have equal probabili- 
ties (= 1/z) for hopping to the left or to the right. The 
RMP, whose position is X(t), is attracted by the gas 
particles and thus has asymmetric hopping probabilities, 
which obey Eq.(2) with ?7(r;t) replaced by p(r;t). 

To determine X(t) we now proceed as follows. Rec- 
ollecting first the results of Refs.9, we anticipate that at 
sufficiently large times the particle density past the RMP 
tends to some constant value. This implies, in virtue of 
Eqs.(l) and (2), that p{X(t)\X(t) ± a) approach limiting 
values p±, which are independent of X(t) and t. Solving 
next the problem for arbitrary fixed p±, we determine 
X(t) and the density profile p(X;t), A being the distance 



from the edge, A = X(t) — X. Finally, inserting p(X;t) 
into Eqs.(l) and (2), we find the closure equation, which 
determines p± (and hence, A) self-consistently as func- 
tions of keT and of the LLI parameters. 

The dynamics of a biased RMP in a ID hard-core gas, 
placed in the just described "shock" configuration, has 
been studied in Ref.ll. It was shown that X{t) obeys: 

X{t) = A v^t, (4) 
where the prefactor A is defined for p_ > p + by 
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$(x) being the probability integral. In the special case 
p_ = p + , A is no longer constant and grows with time as 
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Next, it was found that the density profile p(X; t) past the 
RMP has a characteristic S'-like shape; p(X: t) is almost 
constant (and different from p) in a region of size ~ X(t), 
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while for greater A, A ^> X(t)/A 2 , it approaches p expo- 
nentially fast. Such a form of p(A; t) stems from the fact 
that vacancies propagate into the gas-phase only diffu- 
sively and thus homogenize the density distribution past 
the RMP only at scales of order of X(t). Note, how- 
ever, that the total number of particles is conserved, i.e., 
UniL^ooL^ 1 J L dX p(X; t) = p. 

Turning now back to the 2D problem under study, we 
recall that p± are not arbitrary parameters, but their 
values are determined by the density distribution past 
the edge. Inserting Eq.(7) into the Eqs.(l) and (2), we 
find the self-consistent closure equation for p+ /p- : 

P+ ( a s n P+\ U o( T ) S /q \ 

— = exp(- (3 a 7 e ); 7 e = (1 ) — , (8) 

p_ p_ la 
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where r± denote the 2D vectors (X(t) ± a, Y). 

Therefore, we find that the mean displacement of the 
monolayer edge obeys Eq.(4), in which A is related to 
Uo(T), p and (3 through Eqs.(5) and (8). Below we 
present some analytical estimates of A(e), where e = 
j3 Uo(T) 8/2 is a critical dimensionless parameter. 

We find that depending on the value of the parameter 
e four different regimes can be observed: 
(1) When e € [0;1] the only solution of Eq.(8) is 
p+/p- = 1, which means that here 7 e = and the 
monolayer behaves as an ideal surface gas. In this regime 
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X(t) ~ yjt ln(t) and the density p(A;t) changes rapidly 
with the distance A from the edge being equal to p for 
A — > oo and to zero for A = 0. 

(2) When e €]l;e c [, where e c — —ln(l — p)/p, the pref- 
actor A is positive and finite. Here, X(t) ~ and 
the monolayer also wets the solid. The edge tension 
7 e > and vanishes as j e ~ (T5 — T) when T — > Tb; 
Tb is thus the critical temperature of the surface gas- 
liquid transition, which is defined implicitly by equation 
Tb = Uo(Tb) S/2k B , (e = 1). In this regime A diverges 
when £ — ► +1, A « y/hi(p/ (e — 1)), and vanishes when 
£ -> e c , A « (1 - p) (e c - e)/(1 - (1 - p)e c ). The density 
p(A; t) changes smoothly from the unperturbed value p 
to the value at the edge exp{—(3a^ e ), which is close to p. 

(3) At e = e c the prefactor A is exactly equal to 
zero and the monolayer partially wets the substrate. 
Hence, we denote e = e c as the point of the wet- 
ting/dewetting transition for the monolayer. The corre- 
sponding critical temperature is determined by T w / dw = 
Uo(T w /(j w ) 5/(2 ke £ c ) and depends on the coverage p. 
We note that T w / dw and Tb are simply related to each 
other. When Uo(T) is independent of T (say, for London- 
type LLI) one has T w / dw = T b /e c . Actually, the inset 
in Fig. 2 displays the p-dependence of the ratio T w / dw /Tb 
(= l/£c) for this case. For the Keesom-type interactions, 
when U (T) ~ 1/T, T w/dw = T b /^. We finally remark 
that the relation between s and e c distinguishes whether 
it is favorable, at given physical conditions, to have a 
monolayer with coverage p on the solid surface or not. 
Consequently, knowing e c and the density distribution in 
a sessile drop with the respect to the height above the 
substrate, one can predict the number of superimposed 
layers in the "terraced wetting" regime. 

(4) For e > e c the prefactor A < and the monolayer 
dewets from the substrate. Here, a jammed region (where 
p(X;t) > p, Eq.(7)) of size ~ X(t) appears, which im- 
pedes the motion of the retracting edge; the e-dependence 
of A is thus very weak, being strongly limited by the dif- 
fusive squeezing out of "voids" at progressively larger and 
larger scales. In fact, for sufficiently large e one may ex- 
pect that the dewetting process will be accelerated by 
thickening of the monolayer, as it is suggested in Ref.6. 
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Fig. 1. Initial configuration of a monolayer on top of solid surface, 
landscape created by the solid atoms. 



Wavy lines depict the potential energy 
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Fig. 2. Numerical solution of Eqs.(5) and (8). Solid lines from top to bottom show the dependence of A 
on dimensionless parameter e for p = 0.9, 0.8, 0.7 and 0.5 respectively. The line t = 1 separates the "surface 
gas" and liquid-like phases. The inset displays the p-dependence of e c , (A(e c ) = 0). 
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